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Abstract 

The Gupta-Bleuler quantization method of QED can be gener- 
alized to canonically quantized constrained systems with quantum 
second-class constraints. Such constraints may originate either from 
the second-class constraints, already presented in the classical descrip- 
tion of the theory, or may have their sources in quantum effects, in 
which case the theory is called anomalous. In this paper, I present 
a detailed description of how the Gupta-Bleuler ideas can be imple- 
mented in these cases and I argue that there are in principle no in- 
consistencies in quantum anomalous theories. Having quantized the 
anomalous theories canonically, I derive the path integral formulation 
of such theories and show that some new terms are necessarily present 
in this formulation. As an example, I show how the chiral Schwinger 
model can be quantized in the original fermionic formulation with no 
reference to the bosonized version used in the literature so far. 



*e-mail address: tl8@ nikhef.nl 
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1 Introduction 



The quantum anomaly is the breakdown of classical symmetries of a sys- 
tem caused by quantum effects. Even if the anomalies in global symmetries 
provide us with understanding of an important class of physical effects like 
TT ^ 27 decay, the anomalies of local (gauge) symmetries of a theory have 
been for a long time regarded as an incurable disease. It was often argued 
that such anomalies make it impossible to use the Ward identities, which 
may result in the lack of the proof of renormalizibility, Lorentz invariance, 
and unitarity. For this reason, one of the most often invoked requirements 
for a consistent quantum theory is that the local symmetries of the system 
must be free from anomalies. This requirement has been successfully applied 
e.g., to restrict the representation content under the gauge group of fermions 
in the standard model and to derive the critical dimension in string theory. 

Only few years ago, the applicability of the standard quantization proce- 
dures to anomalous theories has been put in question 0]. It was observed 
that one should not use information concerning the classical symmetry di- 
rectly in the quantum context. The very reason for that are the anomalies 
themselves - the symmetry group of the quantum theory may be different 
from the corresponding classical ones, which is not very much surprising af- 
ter all. Therefore, one may conclude that the problems do not reside in the 
anomalous theories themselves but rather in an inappropriate way of quan- 
tizing them. The problem will become even more severe if one realizes that 
the anomaly may also provide an enlargement of the symmetry group of 
quantum theory^, in which case all standard techniques of quantization (like 
the Dirac approach or standard path integral method) cannot be directly 
applied. 

It seems clear that the only way out of these problems is to discuss the 
gauge symmetries of a theory not prior but after quantization. Since in the 
language of canonical quantization anomalies exhibit themselves as the non 
- closure of the commutator algebra of classical symmetry generators, this 
would amount in developing the quantization scheme such that one can deal 
with a general system of first- and second-class constraints on the quantum 
level. Actually such a scheme is well known for almost 40 years and has its 
roots in the Gupta-Bleuler approach to the quantization of quantum electro- 

^This effect, called antianomaly is briefly discussed in the last section of this paper 
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dynamics 0: 

If one wants to implement the covariant gauge-fixing condition of QED 
dfj^A^^ = directly on states of quantum theory, one encounters an immedi- 
ate difficulty realizing that the negative and positive frequency parts of this 
condition, d^j^A^^^^ = and d^A^~^^ = do not commute and therefore there 
are no nontrivial states satisfying the condition 

d^Af'lphys >= 0. (1.1) 

The simple solution to this problem is to realize that in quantum theory 
the only physically meaningful objects are the matrix elements of quantum 
operators and therefore, instead of ( [1.1|) , one should use 

< physld^Af'lphys' >= 0, (1.2) 

to define gauge-fixed physical state. The idea of Gupta and Bleuler was to 
consider a particular necessary condition for (|1.2| ) given by 

a^^W^lp/iys >= 0. (1.3) 

It was proved that condition ( |1.3| ) describes the physical states of QED cor- 
rectly. 

This method can be easily generalized to the canonical Hamiltonian sys- 
tems with quantum second-class constraints. To this end, one splits the 
system of real (Hermitean) constraint operators Q whose algebra does not 
close 

[Qa,Qb]=ujab (1.4) 

into holomorphic and antiholomorphic parts Qj and Q\, both of whose have 
separately closed algebras and defines physical states as to satisfy 

Gilphys >= 0. (1.5) 

It should be stressed that in this way we effectively treat the holomorphic 
part of the second-class constraints as generators of a gauge symmetry, which 
leaves the physical states invariant. 

In this paper, I describe the Gupta-Bleuler procedure applied to quanti- 
zation of theories with second-class constraints and/or anomalies in general 
setting. After reviewing the standard Dirac procedure in Section 2, in two 



38 



subsequent sections I develop the general Gupta-Bleuler approach. In Sec- 
tion 5, I discuss a possible formulation of the Gupta-Bleuler quantization 
in terms of the path integral technique. In Section 6, I discuss the Gupta- 
Bleuler quantization of the simple anomalous theory, the chiral Schwinger 
model, which in complete agreement with earlier works is shown to be finite 
and soluble and is equivalent to the quantum theory of a massive and mass- 
less scalar fields. Finally, Section 7 is devoted to discussion of the so - called 
antianomaly effect and Section 8 contains conclusions and some remarks. 

2 Review of Dirac theory of quantization of 
constrained systems 

In this section, I give a short review of Dirac theory of constrained systems [^]. 
There are many textbooks in which the extensive discussion of this subject 
can be found: [0, |0| (see also review papers 0, [0); in what follows, I 



will mainly refer to the new excellent book by Henneaux and Teitelboim JlT 
covering also the most recent development on the field. Therefore, below I 
will concentrate only on the aspects of the problem which will be relevant 
for the theory developed later. 

Let the system of interest be described by the Lagrangian L, which de- 
pends on positions qi and velocities g/, where the index / can be discrete 
(mechanical systems) or continuous (field theory): 

L = L{qj,qj). (2.1) 

One defines momenta in the standard way, 

Pi = ^{q,q). (2.2) 
dqi 

It may happen that the Hessian matrix Wr j = ^ = does not have a 

J i-i- J J dqj dqiqj 



maximal rank and therefore not all of the equations ( p.2|) can be solved for 
the velocities. Instead, we have some number of relations 



QA=PA-'ipA{qi,Pi) = 0, (2.3) 
where the maximal number of velocities have been expressed in terms of 



momenta using (|2.2|). These relations are called (primary) constraints, since 
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they restrict the possible dynamics of the system to the surface in the full 
phase space of the system defined by eq. ( p.3|) . 

Given the symplectic structure on the phase space of the problem (i.e., 
the Poisson bracket {-k, -k}), we can calculate the bracket of constraints Qa to 
be in general 

{Ga,Gb} = JaEOc + ooab, (2.4) 

where f^^ and ujab may be phase space point-dependent, and uab is sup- 
posed to be non - vanishing identically on the constraints surface ( p. 31) . We 
assume that uab is such that the Jacobi identity for an bracket algebra ( |2.4| ) 
holds. Among the constraints, there are such that their bracket algebra 
closes, 

{g^,gp} = fjg„ (2.5) 

and which have a vanishing bracket with all other constraints on the con- 
straints surface [Q = 0). In Dirac terminology such constraints are called 
(primary) first-class. All other constraints are called second-class and will be 
denoted by Qi. 

From the Lagrangian of a theory we can construct the canonical Hamil- 
tonian by performing Legendre transformation 

Ho{p, q) = Piqi - L{q, g), (2.6) 

where on the right hand side we replaced all velocities by positions and mo- 
menta wherever possible or use the relations ( |2.3D otherwise. One can easily 
check explicitly that the canonical Hamiltonian defined above is velocity- 
independent. 

Since the dynamics is supposed to take place on the constraints surface, 
we can add to the canonical Hamiltonian a general linear combination of 
constraints and form in this way the so-called primary Hamiltonian 

H = Ho + X^g^ + v'gi, (2.7) 

which is a generator of dynamics. The important consistency requirements is 
that the time evolution takes place on the constraints surface, which means 
that the constraints are supposed to be conserved in time, to wit 

^ = {GA,H} = V^gA. (2.8) 
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It may happen however that the set of primary constraints does not satisfy 
the above equation. Except for uninteresting case when this relation be- 
comes inconsistent (of the form = 1), it is possible that some new relations 
between phase space variables emerge. These new relations are called sec- 
ondary constraints and should be added to the set of primary ones|. Then 
the procedure of splitting the constraints into first and second-class and look- 
ing for new ones should be repeated along the same line of reasoning, until 
all constraints are found, i.e., until eqn. ( p. 81) is satisfied identically. 

The important observation is that the first-class constraints are in one-to- 
one correspondence with the (local, i.e., with time - dependent parameter) 
gauge symmetries of the theory and ( p.5|) is in fact a representation of the 
gauge group. The same is true after quantization: actually the whole idea 
behind the Dirac quantization was to ensure the gauge invariance of quan- 
tum theory. In the Dirac approach, the gauge invariance is introduced by 
demanding that the operators corresponding to the gauge generators annihi- 
late the wave function. The anomaly arises when these operators fail to have 
closed commutators, which means that the classical symmetries are broken 



on quantum level. This important facts are discussed in details in ||Tl . 

Given a full set of constraints Qa, Gi, we can turn to the quantization 
problem. Before doing that, we introduce a concept of Dirac bracket which 
replaces the standard Poisson bracket in the Dirac quantization procedure. 

The second-class constraints satisfy a simple Poisson bracket algebra on 
the constraint surface 

{Gi,Gj}\g=o = uJij, (2.10) 

where Uij is invertible. Then one defines the Dirac bracket of any two func- 
tions on a constraint surface to be 

{A,Br = {A,B} - {A,g,}ico-')Hg„B}. (2.11) 



^It is clear that there is at niost as many secondary constraints as the first-class primary 
ones (keep in mind however that the first-class primary constraints may become second- 
class in the full theory i.e., when all secondary, tertiary, etc. constraints are found as 
described below). Indeed, if Qi is second-class then 

= {H, g,} = {Ho, g,} + VjUj,, (2.9) 

where the last equality holds up to some combination of constraints. This equation is 
always solvable, since w is invertible by definition. 
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This bracket can be understood as a projection of the Poisson bracket 
down to the second-class constraint surface. Thus, we can consistently re- 
place all the Poisson brackets by Dirac brackets and solve the second-class 
constraints, expressing all Pi by ipi using ( |2.3| ). The Dirac bracket provide 
us with a new symplectic structure on the phase space of the problem, such 
that the second-class constraints are absent. The quantization is achieved 
by replacing the Dirac brackets by commutators in the standard way and 
finding the representation of the canonical commutational relations on some 
Hilbert space. 

As for the first-class constraints, the Dirac proposal consists of implement- 
ing them directly on the quantum level. After quantization, the first-class 
constraints become Hermitean operators satisfying the commutator algebra 
(provided there are no anomalies) 

[Ga,G^]=tf^fs'g,, (2.12) 

which is in direct correspondence with the Poisson (Dirac) bracket algebra 



(|2.5| ). In the formula above f^^^ are, in general, operator - valued and, 
of course, one may encounter the notorious operator ordering problem in 
passing from ( p.5| ) to ( p.l2| ). 

Given the operators corresponding to the first-class constraints, one im- 
plements them by demanding that the quantum dynamics takes place on 
the subspace of an original Hilbert space, called the physical Hilbert space, 
consisting of all vectors (wave functions), which are annihilated by the con- 
straints, to wit 

g^lphys >= 0. (2.13) 

As it was said before, since the first-class constraints correspond to the gauge 
symmetries of the theory, this condition means simply that the physical wave 
functions are gauge invariant. 

It follows from the construction discussed above that since the Hamil- 
tonian operator constructed from the primary Hamiltonian commutes with 
constraint operators on physical states, the time evolution maps physical 
Hilbert space into itself. Analogously, we define physical observables to be 
Hermitean operators commuting with constraints. In a final point of con- 
struction we define an inner product on the physical subspace to be the one 
induced from the inner product of an initial Hilbert space of the problem. 
One should keep in mind that very often the initial space of states does not 
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have a positive definite inner product (so, strictly speaking it is incorrect 
to call this space a Hilbert space) and only in correctly constructed theory 
there are no negative and zero norm states in physical sector (if there are 
negative norm states in the physical sector one may try to construct a second 
quantized theory in the standard way). 

Let us conclude this section by discussing the Dirac quantization of a mas- 
sive scalar field theory on x in light cone coordinates^. The Lagrangian 
of the model reads 

L(x+) = — r"(9+<l>9_$ - -m2$2). (2.14) 
27r Jo 2 

Expanding into Fourier modes in x~ 

<l> = 5]$„e^^^", $-n = $L (2.15) 

n 

we have 

H^^) = E H^^nK + m<i>|,$„ - m2$„$t ) - Im^^l (2.16) 

where the overdot denotes differentiation with respect to the light cone time 
. One can readily find the momenta 

: -2n< , (2.17) 



9$„ 

= 7^ = '^^n, (2.18) 

Po = 0. (2.19) 

We see therefore that there are constraints in the theory. The first two of 
them are clearly second-class. The canonical Hamiltonian reads 

i7=^m2$„$t (2.20) 

n>0 ^ 



■^This example will be relevant to what follows. It happens that it describes the effective 
theory of 2 - dimensional chiral QED (Schwinger model) - see Section 6 
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We see that there is a secondary constraints $o = 0, which together with 
Po = form a second class system. We can get rid of them by using the 
Dirac bracket and erase them from the all relevant formulae. Since the non - 
zero 'momentum' constraints (labelled by non - zero index n) do not depend 
on Pq, we need only to modify the Hamiltonian, to wit 

H^J2^'^n^l (2-21) 

n>0 

Let us turn to the remaining constraints 

Gn = P^-^n$„, (2.22) 
Ql = Pn + tn^l (2.23) 

which satisfy the Poisson bracket algebra 

{Gn,GL} = -'2inSn,m. (2.24) 
Now we can write down the Dirac bracket 

{A, BY = {A, B} + ^{A, GnHGl B} - ^{A, GlHGn, B} (2.25) 

use this bracket and replace by $ , P everywhere using the constraint 

equations. The fundamental bracket is now 

{Pn,^mr = -lSn,rn (2.26) 

and the commutator (after rescaling $ — > |$) 

[Pn,<^m] = -tSn,m- (2-27) 

Substituting the constraint (2.30) into Hamiltonian we obtain 

H=Y.m^^^nPn, (2.28) 

n>0 

up to the ordering ambiguity. In the coordinate representation 

^>» = -a|;. (2.29) 
44 



and we have the Schrodinger equation 

^ ^'^^-tI-"^^*") = ^^*(*") (2-30) 

n>0 

Before solving this equation let us discuss the form of the Hilbert space of 
wave functions. In this space the constraint equations must be satisfied on 
the operator level i.e.J^ 

Pt = (2.31) 

where we already rescaled the operators P„ '^^n, \/n^n- This 

rescaling does not alter the canonical commutational relations and the form 
of the Hamiltonian operator. Let us assume that the Hilbert space of the 
problem is a space of polynomials in $„ variables. This assumption follows 
from the fact that the solutions of ( p.30|) are polynomials. For the basic 



elements of such space, the monomials we define the inner product to be 

< a$^, b^i >= a*6(/5!)5„,^5°'^ (2.32) 

It is easy to see that this inner product is consistent with hermicity relations 
(OTp. 



Now, equation (|2.30|) is solvable in terms of linear combination of products 
of finite number of monomials 

^ = []<|.°". (2.33) 

The energy of such solution is 

2 

, 77? 

£^ = J2 (2-34) 

Therefore, we have solved the quantum theory of a scalar field using Dirac 
method. Note that the second-class constraints played an important role even 
after Dirac bracket has been introduced and replaced by the commutator. 
Indeed they provide us with a hint as to how define the Hilbert space of the 
problem and make the Hamiltonian operator Hermitean. One should also 
observe that the Schrodinger equation ( |2.3U| ) does not have any non - trivial 
solution in the space of square integrable functions^]. 

''This shows that the second-class constraints are not totally 'forgotten' by the quantum 
theory even if the phase space of the problem is reduced to the constraint surface with a 
symplectic structure given by the Dirac bracket. 

^It is possible to solve the theory in the Hilbert space being a space of square integrable 
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3 Gupta-Bleuler constraints 



Even though the Dirac procedure of quantization of constraint systems seems 
simple in principle, one may encounter many problems in applications to 
particular models. First of all, the required split of constraints into first 
and second-class may be problematic if one wants to make it covariant with 
respect to some global symmetries of the theory like a global target space 
super Poincare symmetry in superparticle and superstrings models. Secondly, 
it may happen that the algebra of the commutators of constraints 
acquires additional central charge terms, as compared to the Poisson (Dirac) 
bracket algebra ( ^.51 ). This purely quantum effect is in fact the exhibition of 
anomalies discovered originally in the diagrammatic context in late sixties. 
The appearance of anomalous terms in (|2.12 ) makes it impossible to im- 



plement the Dirac procedure of recognizing the physical subspace of the 
Hilbert space as the vector space of states annihilated by constraints. In- 
deed, if the system of operators Qa corresponding to first-class constraints 
satisfy the anomalous commutator algebra^ 

[Ga, Gp] = ifafsG-y + UJa(3, (3.1) 

instead of ( |2.12|) , and we want to follow the Dirac procedure to define the 
physical states as the ones satisfying ( |2.13D , we find immediately that 



^^a/slphys >= =^ \phys >= 0, (3.2) 

since Uajs does not have non - trivial zero modes by assumption. There- 
fore, the Dirac procedure cannot be consistently apphed to quantization of 
anomalous theories and requires modifications. In this sections, I discuss a 
proposal of replacing the Dirac procedure by a new one having its roots in 
Gupta-Bleuler quantization method in quantum electrodynamics. 

There are clearly many methods of imposing constraints in quantum the- 
ory - one needs to restrict a Hilbert space by some reasonable analogue of 

functions. In order to do that, one needs to change the representation of commutational 
relations (2.34) and include an exponential dumping factor in the wave function. Then 
the result is identical to the one presented in Section 4. 

^Strictly speaking, for anomalous theories, the commutator below is defined only with 
respect to some pre - chosen set of states (i.e., one should use the 'sandwiched' form of 
the algebra). Since we always (implicitly) define the Hilbert space as the first step of 
quantization procedure, in what follows we will ignore this subtle point. 
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the classical condition of vanishing of some combination of phase-space vari- 
ables. Among these possibilities, the Dirac requirement that the constraints 
annihilate the physical states, ( p^.l3D is probably the most restrictive. As we 



have shown above, it does not work for the system of mixed constraints, and 
we must get rid of the second-class ones on the classical level by making use 
of the Dirac brackets. The other, physically more appealing possibility is 
to define physical states as the maximal subset of all states such that the 
matrix elements of all operators corresponding to constraints vanish between 
any two of them, to wit[| 

< phys\gA\phys' >= 0. (3.3) 
Clearly, in this case the constraints operators Qa can form the algebra 

[Qa, GB]=tfAEGc + ^AB, (3.4) 

with the non - vanishing central charge Uab which corresponds either to 
the classical second-class constraints, or has its source in anomalies, without 
breaking the consistency. However, condition ( |3.3| ) is not very restrictive and 
one has problems with interpretation of physical states defined in this way 
(cf. discussion in [p!2|). 

In this paper, we choose a middle of the road method. Suppose that 
we can split the full algebra of Hermitean constraints Qi into two subsets 
consisting of the complex constraints Qi and their Hermitean conjugate Qj 
such that the algebras of holomorphic and antiholomorphic constraints close 

[gI,gJ] = ^fI/gK, (3.5) 



=-^G]<nJ, (3.6) 



Tt IK 



gj,g\ =2Zjj. (3.7) 



Then, one defines the physical states as a subset of the solutions of (|3.3| ) 
defined as 



^It is implicitly assumed that the Hilbert space of the problem does not include any 
zero-norm states \ipo >G TCq, || |V'o > ||^ = 0- If such states were present initially, they 
should divided out: H ^ H/Ho- 
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Gilphys >= 0. (3.8) 



Clearly, if (|3.8| ) is satisfied, the condition (|373|) is satisfied as well - indeed it 
follows from ( |3.8| ) that the matrix elements between physical states of both 
Qj and Qj are equal to zero. 

The important question arises, as to whether the Gupta-Bleuler algebra of 
constraints ( |3.5| - |3.7| ) is equivalent to the original algebra commutator ( |3.4| ). 
Before giving the detailed proof of the equivalence of these two presentations 
let us describe how (|375| - |3^ ) can be constructed from ( |3.4| )p|. 

The procedure consists of two major steps: the disentanglement of the 
first-class subalgebra and the polarization of the remaining subalgebra of 
(|3.4| ). Let us stress already at the beginning that both steps may in general 
cause breakdown of some global symmetries present in the theory and that 
some of resulting operators may be non - local. 

In the first step, we construct the maximal set of vectors labelled by the 
index a, such that v are built out of Hermitean operators and they commute 
both with Qa and uab- Suppose the following condition is satisfied: 

V^COAB = 0. (3.9) 

Then it is easy to see that Qa = v^Ga are first-class i.e., they form the closed 
algebra with all constraints. Indeed 

Va^BA + terms proportional to Qa, (3.10) 

and from ( |3.9| ) the first term vanishes. Observe that by construction Qa are 
Hermitean. 

Now we split the whole set of constraints (Qa) into two subalgebras, one 
formed by Qa and defined above and the second consisting of remaining 
operators, which we will call Qi, which satisfy 

[Q^,QJ]=^f^,''GA + UJ^J, (3.11) 

where now the operator- valued matrix uJij has a maximal rank. The polariza- 
tion of the second-class constraints Qi amounts in constructing an operator- 



^This question has been addressed before in ml 
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valued matrix Jij (a complex structure), which commutes with Qi and satisfies 

JijJjk = —Sik, (3-12) 

Jij^jk ~ Jkj^ji = 0. (3.13) 

It is easy to see that once these conditions hold, the complex constraints 
Gi = Qi + iJijQj and Q\ = Qi — iJijQj satisfy the algebra of the form ( ^75| - 



|3.7| ). In the final step we form a full set of generators of the Gupta-Bleuler 
constraints algebra by taking Qi = {Qa, Qi) and Q\ = {Qa, Qj)- 

One should observe that the constraint algebra can be transformed into 
even simpler algebra, namely 

[G«, G^] = [G„, Gi] = [Gi, G,] = 0, (3.14) 

[Gi,G]]=tSi,, (3.15) 

which I will call 'abelian Gupta-Bleuler form' (here and below I denote con- 
straints in abelian form by G). It is clear from this form of the algebra that 
the holomorphic and antiholomorphic constraints essentially behave like an- 
nihilation and creation operators. Indeed, it follows from ( p. 14 ), ( 3.15 ) that 



there exist a local coordinate system (P,Q;pa,qa;pi,qi) in vicinity of the 
constraint surface such that there exist a symplectomorphism G^ —>■ Pa, 
Gi ^{Pi + iqi) G| \{pi — i(li)- Here P and Q are coordinates 'orthogonal' 
to the constraints surface. 

As we saw above, the problem of bringing an algebra of constraints to the 
(abelian) Gupta-Bleuler form is a problem of solving systems of differential 
equations for v and J. Using Darboux theorem, it is possible to demonstrate 
(cf. [|1^ and the detailed discussion in |jl^) that these equations are always 
solvable locally^. One should stress however that the proof is essentially 
classical in the sense that it make use of notions of phase space, constrained 
surface, and Poisson bracket. Therefore, one should be very careful about 
resolving the ordering problem at the beginning of the procedure and tracking 

^The only exception is a system consisting of odd number of fermionic constraints. 
Subsystems are pretty rare and, besides, one can applied the Gupta - Bleuler theory to any 
subsytem of even number of constrains and solve the remaining ones by the standard Dirac 
method. Moreover, since our main goal to construct the consistent quantum theory with 
anomalies, such problems never arise (systems with odd number of fermionic constraints 
are quantum mechanical not field theoretical). 



49 



any ordering ambiguity which may arise so that the correspondence between 
commutators and Poisson brackets is kept. 

On the other hand, it should be observed that the proof merely guarantees 
solvability of the polarization problem and is not operational in general (in 
the sense that it does not prescribe the step - by - step procedure). This 
procedure is either obvious (as in the case of examples considered in this 
paper) or very difficult (as in the case of anomalous Gauss law algebra of 
chiral QCD, see |15|). In the second case, it is little hope that the canonical 



approach is capable of producing anything more than a perturbation theory; 
but then one can equally well turn to the path integral approach, which, as it 
will be shown below, does not require explicit Gupta - Bleuler polarization, 
but only its existence. 

Therefore, the conclusion is that the Gupta - Bleuler procedure is feasible 
in all cases of interest. 



4 The Gupta-Bleuler Quantization Method 

Let me now describe the Gupta-Bleuler quantization method for theories 
with second-class constraints and/or anomalies. We start with a classical 
Lagrangian L. From it we deduce a phase space with its symplectic struc- 
ture, the set of primary real constraints Q\ constraints and the canonical 
Hamiltonian Hq. Then we quantize: we define a representation of canonical 
commutational relations 

[Pi,qj] = -i5iJ (4.1) 

on some Hilbert space (again, strictly speaking it will be an inner product 
space since some vectors may have negative norm). From the algebra of 
Hermitean constraints, 

we extract the first-class ones, Qai polarize the remaining ones, and obtain 
in this way two sets of holomorphic and antiholomorphic constraints, 
where both of them include the first-class ones. In this way the algebra of 
constraints is cast to the Gupta-Bleuler form (|3^ - p.7|) We define primary 
physical states \phys, > to be annihilated by the primary holomorphic 
constraints, to wit 

g'}\phys,0>=0 (4.3) 
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It is clear that in theories with anomahes it would be inappropriate to look 
for secondary constraints already on classical level. Indeed, as was noted in 
footnote 2, the secondary constraints results from primary first-class ones 
and since the number of them may change due to the quantum effects the 
number and form of secondary constraints may change also. 

We start the construction of secondary constraints with definition of pri- 
mary Hamiltonian operator if, which by assumption generates a time evo- 
lution by virtue of Schrodinger equation, 

H = Ho + u'g^j + g^jK'. (4.4) 

Observe that H is Hermitean between primary physical states. In complete 
analogy with Dirac approach, we assume that the time evolution maps physi- 
cal states into themselves, and therefore the primary Hamiltonian ([4.4| ) must 
commute with holomorphic constraints on physical states. We have therefore 
an equation 

= [Q^,, H]\phys, >= {[g'j, Ho] + Z^v' + G'JiG'j, v'])\phys, >, (4.5) 

which after little algebra can be rewritten as 

{g'jHo + g^g°j^v')\phys, O >= O. (4.6) 

The equations above can be either solved for , or produce new constraints 
which should be included into the constraints set, polarized, if necessary and 
then the procedure of hunting for new constraints should be repeated. Af- 
ter the whole set of constraints of the theory has been revealed, one can 
construct the observables in complete analogous way: take any Hermitean 
operator, add a linear combination of constraints and then check if it com- 
mutes with holomorphic ones on physical states. Having obtained the full 
set of constraints ^7 and g\ we define the physical states to be annihilated 
by the holomorphic ones and the time evolution is subject to Schrodinger 
equation 

{^^-H)\phys,t>=0, (4.7) 

where H is defined in (|4.4| ). This concludes the construction of quantum 
theoryQ. 

^"^Let me stress here again: the idea was to take as few bits of information from the 
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Let us now return to the example of Section 2 and analyse it again from 
the point of view of Gupta-Bleuler method. 

After performing the same steps as before, we see that the system is 
described by the set of holomorphic and antiholomorphic constraints 







(4.8) 


Go 


= Po - irri^^Q, 


(4.9) 




= P„ + m$|,, 


(4.10) 






(4.11) 



and the canonical Hamiltonian 

H^Y.^'^^n^l (4-12) 

n>0 

The Gupta-Blculcr Hamiltonian can be obtained by using the procedure 
described above and reads 

HgB = E(^n< + -^nPn). (4.13) 

^ n>0 ^ 

Let our Hilbert space be a space of square integrable functions ^. Solving 
the holomorphic constraints Qn^ — Go^ — ^ gives 

^^g-^*§g-E„>o"*"*^i($„), (4.14) 

where is a polynomial. Substituting this result to the Schrodinger equa- 
tion 

\m' + -*nPn))^ = (4.15) 

^ n>0 ^ 

we obtain, after subtracting the infinite energy due to the ordering change 
in the second term an equation for 



iiv? 



2n 



$„P„^ = (4.16) 



classical theory as possible the primary constraints and the primary Hamiltonian and 
perform all subsequent steps on quantum level. Of course, one can try another approach, 
namely, to reveal all classical constraints and canonical Hamiltonian on the classical level, 
then quantize and try to make the resulting theory consistent as described above. There 
is little hope that these two procedures give the same result in general. 
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which is nothing but the equation ( |2.3CI| ) discussed in the Section 2. Therefore 
we conclude that the Gupta-Bleuler approach gives rise to the same energy- 
spectrum and Hilbert space as the standard Dirac procedure[^. One should 
observe that the reason for using holomorphic constraints is that there are 
no solutions to the antiholomorphic constraints equations in the space of 
integrable functions. Our construction above can be readily applied to the 
massless case, = 0. In this case the Hamiltonian is zero and the theory is 
described by the holomorphic constraints P„ — in^n and the single first-class 
constraint Pq. The wave function of the massless scalar field is therefore 

^massless = S">0 V]/, ($„) . (4.17) 

To show how the Gupta - Bleuler technique work in the real field theo- 
retical framework, let us consider yet another example, namely the massive 
four-dimensional QED. The Lagrangian of the theory reads 

L = -^F.^F^-' + Kn'A.A^ = 

= -\{F^,f + \{d,A, - d,A,f + ^m'Al - \m'Al (4.18) 
By using standard procedure, we easily find two sets of constraints 

Qq{x) = diPi{x) + mMo(x) ^ 0, 
Poix) ^ 

and the Hamiltonian 

H = l d'x (^^P^ + \m-\d.P,f + \m'Ai + 1^ + XP, + pQo) • 

We assume the following form of canonical commutational relations 

[P^{x), Ay{y)] = iri^J^{x - y), r/^^, = (1, -1, -1, -1). 

The constraints Qo and Pq are clearly second class. Since Qo does not 
commute with Ai, we make the replacement Ai —>■ Ai = Ai + m~^(9jPo- It is 

^^The only difference in the form of the Hilbert space, the e~~ » factor, is irrelevant, 
as the only observable which does not annihilate it is a constant in $0: sector. 
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easy to see that the change of Hamiltonian can be absorbed into redefinition 
of A so that we have 

H = j (fx {^Pf + ]^m-\d,P,f + ]^m^Al + ^4 + APo + pQi 



with 



[Pi{x),Aj{y)] = -SijS\x-y), 

[Po{x),Qo{y)] = im'6\x-y). 

Dirac quantization of the theory is strightforward. One just use the Dirac 
bracket and forgets about Pq and Qo whatsoever; it remains to look for energy 
eigenstates of the Hamiltonian 

Interestingly enough, the theory is fully solvable. To see that we first 
remove the ^F^j term from the Hamiltonian by replacing 

^^J^(A) = e^f'^'''''''-^^^^-^^^/j{A). (4.19) 
The wave function ip{A) satisfies 

J d'x (^^P^ + ^m-\d,P,y + l^'Al) ij{A) = £^ij{A). 

It remains only to decompose the fields and momenta into longitudinal and 
transversal parts 

Ai = Aj^ ■j==.A , QiAj^ 0, 
V A 



[P^{x),Aj{y)] = -z6,6'{x-y), 
[P\x),A'^{y)]=^VAS%x-y) 

and realize that the Hamiltonian is a sum of independent three oscillators 
(the term in exponent in ( |4.19|) depends on A"^ only): 



H = Jd'x {(IpP + \m-Ar) + {m-^\pl' + \m^Af': 
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Therefore the wave function resulting from the Dirac quantization is the 
wave function of three independent oscillators (Hermite polynomial times 
the exponential dumping factor)^ multiplied by the term ( |4.19| ). We will 
not discuss here the problem of construction of the inner product. 

The Gupta - Bleuler procedure is straightforward. One takes the holo- 
morphic constraint to be Pq — iQq and realizes that the wave function as 

Jgb(Qo,A) = e-^/«o^D(^), 

where \E'_d is the wave function described above. We see therefore that the 
Gupta - Bleuler procedure merely re-introduces the variables which have 
been cut out classically in the Dirac procedure and the addition part of the 
wave function takes the form of the oscillator vacuum. In view of the fact 
discussed in the previous section that Gupta - Bleuler constraints can be 
identified with annihilation and creation operators this result is very natural 
and indeed comprise the essence of the Gupta - Bleuler procedure. 

To finish this section let me draw the reader's attention to the following 
important fact. Contrary to the Dirac procedure where the physical states are 
not normalizable (as the wave function does not depend on some variables) in 
the Gupta - Bleuler procedure, as it was discussed above, the physical wave 
function contains exponential dumping factor and the physical states are 
perfectly normalizable. This circumstances make it possible to interpret the 
inner product of physical states directly in terms of physical probabilities. 
Therefore, one avoids the neccesity of introducing projectors on physical 
states, which is one of the major technical problems of Dirac quantization[^ 

5 Path integral 

The path integral technique (for review see e.g. [0) is nowadays the most 
popular way of quantization. The simple reason behind it is that it serves 
the most economic way of obtaining Feynman rules (i.e., the description of 

^^The argument for the longitudinal part is A^A^^A^. 

^•^One could think about overcome this problem as follows. Given Dirac (first class) 
constraints, add gauge fixings, polarize the system and follow the Gupta - Bleuler pro- 
cedure. This proposal is certainly worth further investigation, however in many theories 
finding the correct gauge fixings is a very difficult problem itself. 
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perturbative sector of quantum theory) directly from classical Lagrangian - 
it amounts in formal manipulations of the vacuum-to- vacuum amplitudes 



Z[J] = J V^iexp{i{S,i + J0)}, (5.1) 

where Sd is a classical action, and the sources J are coupled to all fields 0, and 
P/i is an appropriate measure. Complementary to canonical quantization in 
which the time evolution of wave function is unitary by construction but 
some global symmetries of the theory are not obviously preserved, in path 
integral all symmetries are manifest^ and one needs to check unitarity of the 
S-matrix|^. 

It is not easy if not impossible to derive the expression ( |5.1| ) directly 
from the canonical analysis in general case. What one can do more or less 
rigorously is to express the matrix element of the evolution operator in terms 
of the path integral over the phase space for unconstrained system, to wit 

<g',t'|g,t> = < g'|e^(*-*')^|g > 

= / \{dq7:-exp{i dt{pq - H{p,q))}. (5.2) 

J path ^ ATI Jt 

Then one can try to integrate over momenta to arrive at the formula ( |5.1| ). 
This is however possible only if the integral is Gaussian. 

Let us pause here for a while to see how formula ( p. 2] ) can be derived and, 
what is more important, to understand what is the meaning of the symbols 
that appear in it|^. 

Let us consider the probability amplitude < (f,t'\q,t > (we use vectors 
over the symbols to remind the multicomponent nature of \q >, \p >). Due 
to the principle of unitarity, the time interval {t' — t) can be split into + 1 
subintervals, each of the length e and the first equation in (|5.2| ) can be rewrit- 



^"^As observed by Fujikawa [||, even though Sd is invariant with respect to aU relevant 
symmetries of the problem by construction, the measure, Vfi may change by the non - 
trivial Jacobian. This is the way how the anomalies show-up in the path integral language. 

^■^ Obviously, this holds true for Lagrangian path integrals. The Hamiltonian path in- 
tegral which I will consider below has the same virtues as the canonical quantization 
procedure. 

-"^^This subject is extensively and deeply discussed in the book by Berezin and Shubin 
ilql and the reader can find there lots of details which we omit here. 
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ten asQ 

< q,t'%t >= 

= J iUd^i) < (fle-'^'^WN >< qN\e-''''\qN-i > ■ < gi|e-^^^|g >(5.3) 

This formula is still exact, but the presence of operators in the exponents 
makes it not very useful. The idea is to replace the matrix elements of 
operators by some other object, which is easy to calculate. To make this 
crucial step, we need to introduce a notion of the classical symbol of operator. 
First of all, we must resolve the ordering problem by defining an expansion of 
any operator A into momentum and position operators p, q, respectively (for 
a while we introduce hats here to distinguish quantum- mechanical operators), 
to wit 

Taking the position and momentum eigenstates, \q >, \p >, < g|p >= 
exp(ig ■ p), one can easily find that the matrix element of A equals 

A{p, q) =< ^A\p> e""^-^, 

where, 

~ 5Z •^h-iaji-j/sPji ' ' 'Villi ^ ' ' ' Q " ■ 

Now we are in position to define the path integral, eq. (|5.3| ). We have, 

< Qtle'^'Wj >= 
dp 



where Ti is the symbol of the operator H and 0{e^) terms reflect the difference 
between the symbol of exponent and exponent of the symbol. 

Now, it is customary to take the limit e — > 0, Ne = t — t' and denote the 
result as 

^^-iS(pq-H{p,q))dt ^ I [dp][dq]e-'^^'. (5.4) 

Some important remarks are in order 

""^^We assume that the Hamiltonian H is not expHcitly time-dependent. 
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1. The integration in ( |5.4| ) runs over the set of eigenvalues of momentum 
and position operators in their eigenstates. This shows again that the 
path integral makes use of the semiclassical and not classical objects. 

2. The q in the final formula is a shorthand notation for the expression 
lime^o ■ But, as it is well-known, the major contribution to the 
path integral comes from the trajectories, which are not differentiable. 
One should also carefully check if the 0{e^) terms do not contribute to 
the path integral. One may also understand this term as reflection of 
inherited ambiguity in definition of the path integral which disappears 
if one sticks to one well - defined prescription of how path integral is 
constructed. This subtle point is discussed in [|l^]. 



3. Even more important is the simple but often overlooked fact that the 
terms in the exponents of path integral are the symbols of operators 
and not the classical objects. These two may well coincide, but there 
are important examples where they do not. Actually, the anomalous 
theories serve as the primary examples of such situation. The point 
is that if the Poisson bracket appears in the path integral exponent, 
one should understand it as the quantum mechanical matrix element 
of the corresponding commutator, between position and momentum 
eigenstates (it is easy to see that instead of momentum and position 
eigenstates one can take any other complete set of eigenstates as, for 
example, the particle number operator eigenstates (Fock basis)). 

Keeping the derivation and comments in mind, let us now turn to dis- 
cussion of the path integral quantization of gauge theories. Below, we will 
discuss only the phase space path integral which, contrary to the Lagrangian 
path integral, is directly linked to the most fundamental canonical quantiza- 
tion technique[^. 

For constrained systems, the theory of phase space path integral was 



developed by Faddeev [|I^ in the case of first class irreducible constraints 
and generalized by Senjanovic to incorporate the second-class ones 0. The 
main idea was simple to postulate a form of the path integral by adding the 



^^The problem with the Lagrangian path integral is that, in general, there is little control 
over theory being quantized, especially when problems like anomalies arise. 
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delta functions of constraints and gauge fixings together with appropriate 
determinants to the measure and integrate over the unrestricted phase space 
of the problem. Then one checks if the resulting expression is the same as one 
obtained from ( p.2|) for the integration over the reduced phase-space of the 
problem (i.e., when the constraints are solved by means of (classical) gauge 
fixings prior to construction of the path integral). The expression ( ^.1| ) can 
be then obtained after integration over momenta as before. 

In the middle seventies, in the remarkable series of papers |19[, |2^, ^1 



see also reviews ||22|, [^) the most general form of phase space path integral 



was obtained by means of extending phase space as to incorporate ghost 
fields, which have had been introduced earlier in Lagrangian formulation of 
path integral by Faddeev and Popov |Q. In this extended phase space the 



local gage symmetries of the problem are replaced by single global symmetry 
called the BRST symmetry. This symmetry was discovered in the middle 70's 
as a symmetry of Faddeev - Popov path integral and it was realized later 
that it plays important role both in classical and quantum theory. Since our 
formulation of path integral is essentially based on the BRST construction, 
let us pause for a moment to review it. 

The volume of this paper does not permit even a brief overview of the 



BRST theory, the reader should again refer to [|TT] for more detailed infor- 
mation. Let us just recall the most basic notions. 

Given a set of (Hermitean) first - class constraint operators Qj^ 

one defines the nilpotent BRST operator[^ 

Q = C-Ga - Ic^c^flpTT,, (5.5) 

[g,g]+ = o. (5.6) 



^^For simplicity I assume that all constraints are bosonic, irreducible, and their algebra 
is not open (/J^ below are constants). The most general case is discussed in [0. 

^"^In quantum mechanical framework there is of course the ordering problem which must 
be resolved. Here and below, I assume that some ordering prescription (in our case normal 
ordering) is always in force. Otherwise, one should add some additional terms, e.g., to 



make (5.5) Hermitean. 
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Since the BRST charge Q is nilpotent, it provides a well-posed coho- 
mology problem. To proceed, we introduce a grading in the space of wave 
functions over extended phase space as follows: We say that the 'ghost num- 
ber' of the ghost c is +1 and of the 'ghost momentum' —1. We assume 
further that the wave functions are (finite) polynomials in ghosts and there- 
fore, they can be decomposed into a sum of terms with definite ghost number. 
Acting on such space, the BRST operator is a cohomology operator, since 
it raises the grading by one. We define the BRST (state) cohomology by 
|-|(") = Ker Q\n/ ImQ\n-i (the kernel of Q acting on the space with ghost 
number n divided by the image of Q acting on the space with ghost number 
n — 1). Then it is obvious that the cohomology space with no ghosts is the 
physical space of Dirac. Indeed, KerQ\Q = fla KerQa and Jm = 0. 
It is assumed that the generator of dynamics has zero ghosts number and 
commutes with the BRST charge, therefore, the cohomology does not change 
during the time evolution and if one assumes that the initial state is physical, 
it will evolve to a physical state. 

It is not necessaryin general that the BRST charge has the form ( |5.5| ). In 
fact, this is the so-called minimal form of this operator. One can add any 
additional sector with trivial cohomologies, for example, the sector containing 
the Lagrange multipliers A" with their momenta vr^ and the corresponding 
ghosts sector a a- The additional charge 

Q = PaTT 

has clearly the zero cohomology being the wave functions independent of A 
and the zero cohomology space of Q + Q"''^'^ is 

sfHqlo (g) {functions independent of A}. It turns out adding that such 
a BRST charge is very convenient in formulation of the path integral. For 



more detailed discussion, see 11 



There is one more important requirement concerning the BRST charge. 
In order to assure automatic decoupling of closed states \ip >= Q\k > from 
physical states, the BRST charge must be Hermitian. Indeed, then 

< phys\ip >=< phys\Q\K >=< K\Q\phys >= 0. 

Moreover, for Hermitean Q, the closed states have zero norm. In other 
words, the hermicity of the BRST charge with respect to some inner prod- 
uct guarantees that the inner product can be consistently restricted to the 
cohomologies. 
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The system of constraints {Ga} is, of course, not unique, as it can be 
replaced by any other system {A^^/j}. One of the important properties of 
the BRST formahsm is that such transformation is generated by a canonical 
transformation in extended phase space and therefore it leaves the measure 
of the path integral, which is a Liouville measure by construction, formally 
invariant. 

The proof of this fact consists of two steps. Since is invertible by 
definition, det 7^ 0. If det A^ < 0, by making the canonical transformation 
of ghosts ^1 ~^ ^-iid all other ghosts invariant, we get the 

BRST charge corresponding to the system of constraints {—Qi, Ga}, {(^ > 1)- 
Therefore, without loss of generality we can assume that det A^ > 0. Now, 
A(^ can be obtained as a composition of infinitesimal transformations 

= 5f + e^. (5.7) 
These transformations are, in turn, canonical transformations generated by 

Let us turn to discussion of dynamics of the theory. The BRST analogue 
of the fact that constraints are preserved in time is that the BRST charge 
is time-independent. Therefore, one looks for extension of the canonical 
Hamiltonian Hq which commutes with Q, to wit 

H = Ho + ..., [H,Q] = 0. (5.8) 

Such a generator of time evolution is called the unitarizing Hamiltonian. It 
is easy to observe that the unitarizing Hamiltonian is not uniquely defined 
by eq. (refunitham) . Indeed, one could replace 

H ^ H +[K,Q] + . (5.9) 

In the formula above, the function K which is assumed to have the ghost 
number —1, is called the gauge fixing fermion. In fact, it can be shown (see 
[pip that different choices of K correspond to different gauge fixings condi- 
tions. 

Let me now summarize all the previous results and present the step-by- 
step procedure of construction of the Hamiltonian path integral which I will 
follow below: 
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1. 



In the first step one constructs an extended Hilbert space built over 
the phase space which includes both physical and ghost variables. 



2. Next, one constructs the BRST operator acting in this Hilbert space 
and being nilpotent and Hermitean. The zero cohomology space of 
this operator is assumed to coincide with the space of physical states 
of Dirac quantization. 

3. Given the BRST operator, one chooses the gauge fixing fermion K and 
constructs the unitarizing Hamiltonian operator (|5.8|) . This operator 
plays arole of the generator of dynamics. 

4. The unitarizing Hamiltonian is used to derive the path integral follow- 
ing the formal procedure explained above. It is important that any 
formal manipulations on the path integral can be made only after the 
symbols of operators are computed, i.e., when the objects under dis- 
cussion are already well defined functions on the classical phase space 
of the problem. 

Let us employ this procedure in the case of Gupta-Bleuler constraints. 
Our first problem is to find an appropriate BRST operator. This step 
has been already performed in the paper |Jl3||. It was shown in this pa- 



per that one can define the so-called 'second order BRST operator' Vt which 
is Hermitean by construction and whose zero cohomology coincides with the 
Gupta-Bleuler physical states. This operator is constructed as follows 

Vt = s^Q + QU + ^Z -2shl3, (5.10) 

where (s^, s) are second level commuting scalar ghosts and (7, (3) is a pair of 
real conjugate Fermi-type second level ghost operators: 

"s,st] = 0, (5.11) 
{7,/5} = l, (5.12) 



Q = c^Qa + c^c^fsA^c + ... (5.13) 
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= c^^G\ + -ywc^^ + (5.14) 

{Q,Q^} = 2Z, (5.15) 



[Z,Q]= Z,Q 



0. 



All details of construction of the inner product in the Hilbert space of the 



problem were presented in . It should be only mentioned at this point that 
the cohomology of Q reduces to a sum of (dual in some sense) cohomologies of 
Q and . The main result of this paper can be summarized in the following 
diagram 

Q^Q I (5.16) 
V* ^ V{ ^ ... ^ 0. 

In this diagram the space V2a is first level ghost free and corresponds to the 
physical space of the Gupta-Bleuler procedure. Moreover, the spaces Vi and 
V* are adjoint with respect to the inner product of the Hilbert space of the 
problem. 

In the next step, we construct the unitarizing Hamiltonian. Let us assume 
for simplicity that the canonical Hamiltonian Ho commutes with both Q and 
Q^. There exists a very natural choice of the gauge fixing fermion K which 
is simple and leads directly to the desired form of the path integral, namely 
K = (3. Indeed, for such choice the unitarizing Hamiltonian reads 

H = Ho + [Q\Q] (5.17) 

Let us observe, that in the formula above Q'^ can be interpreted as a gauge 
fixing fermion for the Gupta-Bleuler BRST operator Q, which is a very 
natural choice indeed. 

At this point we can forget the second level ghost whatsoever. One can 
think that they were integrated out and are no longer present in the con- 
struction. 
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As it was said before, it is convenient to extend the BRST operator by 
adding an additional sector of variables. Thus, let 

Q^Q + PaA^^, Q^-^Q^ + p\\^ (5.18) 

At this point, one must further specify properties of the operators which 
are introduced in definitions of Q and Q^. First of all, we identify the variables 
c^^ with the momenta of the ghosts pA (and accordingly for their conjugate 
variables). Next, we define A"^ to be momenta of A^^: [A"^, A^"^] = —i. Then 
it is easy to see that the zero cohomology of Q is a tensor product of the 
Gupta-Bleuler physical states and the states independent of the Lagrange 
multiplier A^"^ enforcing the constraints Qa- 

With these definitions in hand, we can further simplify formula ( ^.17] ) 
and the form of path integral. However, here one should proceed with care. 
The reason is that there are examples of theories (like strings) which have 
anomalies in the ghost sector. For such theories, some of commutators in 
( |5.17D may have anomalous terms and such terms, if present, could not be 



omitted. However, for simplicity, we will ignore this possible complication 
(One of the applications of the formalism presented here is to formulate the 
path integral for string theory away of critical dimension. This problem will 
be discussed in the separate paper.). Keeping this remark in mind, we can 
proceed and just compute ghost commutators. For the sake of simplicity of 
the following few formulas, I will assume that the constraints Q and Q'^ are 
in the abelian form - the reader can easily convince himself that the final 
formula (^.19|) below remains the same in the general case. 



Thus, we have 

Q = c^Qa + Pa\'^^ 

and the final form of the symbol of the unitarizing hamiltonian to be inserted 
into the formula for the path integral (p.4|) reads 



n = no + c^^'iGl Qb\q c"" + \^^Q\ + \^Qa + ip\pA. (5.19) 

In the formula above we introduced the notation [*,*]q to denote the 
symbol of the commutator. Clearly, for the theories with anomalies this 
symbol does not coincide with the Poisson bracket and contains the additional 
anomalous terms. 
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There is one more important group of terms in the Sd in ( p.4| ), namely 
the kinetic terms. Apart from the kinetic terms for classical phase space 
variables we have (in agreement with our choice of hermicity properties of 
this sector) 

Skin = p\c^ - c^^PA + ^At^A^. (5.20) 

To transform the formulas above to a more familiar form, we can perform a 
couple of (formal) manipulations on the ghost and A sector of the resulting 
path integral. First, we integrate over p and and then we make a (singular) 
scale transformation: 

Q\^^-^Q\^ ct^^ect^, A^^eA^. 

This transformation does not change the measure of the path integral, 
since there is as many A as c variables and they have opposite statistics. 
Now, we can take the limit e — to obtain the final form of the path 
integral 

Sgb = pq-no- c^^'iGl Gb\q - At^e;! - \^Ga. (5.21) 

Z = j[dp][dq][dc]d[c^][d\]d\^]e-'^^^ (5.22) 

As it was mentioned in Section 3, is usually quite hard to polarize con- 
straints and this is one of the major problems of Gupta - Bleuler procedure. 
The remarkable fact concerning the path integral ( ^.22|) is that one does not 
need to do that. In fact, one can start with path integral with holomorphic 
'constraints' and antiholomorphic 'gauge fixings' and go all the way back to 
real constraints. 

As it was discussed in Section 3, one of tha basic properties of the Gupta- 
Bleuler construction is that there exists a one - to - one linear relation 
between real and holomorphic (antiholomorphic) constraints, to wit 

Ga = J^\Gi, G\ = ^AGi- 

Then, we have 

(A^e^i + X^^Ga) = (A^Ajj + At^A^)^;, = yG^. 
As far as the ghost term is concerned, 

AgIGb]qc'' -^c%,G,]qc^. 
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In the formula above, the terms proportional to Qi were included into X^Qi. 
Now it is easy to observe that the Jacobians of transformations A"^, A^"^ — >■ A' 
and c"^, c^^ d are the same but since these fields have opposite statistics, 
they cancel in the measure. Thus, the path integral for anomalous theories 
can be written as 

Z = [ [rf/i]e^^-'+'="'^'^='. (5.23) 



where Sd is the standard classical action and the effect of anomalies reside 
in the additional ghost term. 

It is clear that the formula ( p.23[ ) can be also applied in the case of 



theories with classical second class constraints. In this case this formula 
coincides with the formula derived by Senjanowic ^ and this fact can be 
treated as the first check of the above formalism. 

It is interesting to note that the term Uij above is known in the case of 
four dimensional chiral QCD and one can try to use ( p.23| ) as a departure 



point for discussion of the perturbative quantum theory of anomalous four 
dimensional chiral QCD. 

To finish this chapter, let us observe that the formula (|5.23|) clearly differs 
from the standard path integral formula for gauge systems. In the latter case, 
the path integral reads 

\dfj]e''^^'+'''^^-^^^'^. (5.24) 



The differences are twofold. First, the formula ( |5.24| ) contains twice as many 



ghosts (c, 6) as (|5.23| ). Secondly, according to the standard formalism, in 



( |5.24| ) one introduces the gauge fixings JF for the gauge generators Q. It is 
obvious that such approach leads to inconsistencies since Q are not generators 
of symmetries in the quantum theory if anomalies are present. 

The result ( p.23| ) can be also understood as follows. In the presence of 
anomalies, the gauge symmetries of the system is reduced by half (as it is 
clear from the Gupta-Bleuler picture) and the anomalous constraints start 
playing a role of their own gauge fixings. 



6 The chiral Schwinger model 

In this chapter, I demonstrate how the ideas developed above work in prac- 
tice. There is not hard to find an example for this demonstration - indeed the 
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chiral Schwinger model provide us with the theory which is both simple 
(actually quantum mechanically solvable) and interesting as a toy model for 
four dimensional chiral QED or QCD. 

As an anomalous theory, the chiral Schwinger model has been for long 
time considered to be quantum mechanically inconsistent. Only in 1984 



Jackiw and Rajaraman, ||2^, discovered that this model can be consistently 
quantized and showed that the physical spectrum consists of the massive and 
massless scalars, with mass parameter, m, proportional to the electric charge 
of the original theory, e. However, Jackiw and Rajaraman did not refer 
explicitly to the model in its original fermionic formulation, but rather they 
considered the model with bosonized fermions. This formulation has a virtue 
that there are no anomalies: the anomalous commutator in fermionic formu- 
lation was replaced after bosonization by the classical second-class constraint. 
The serious disadvantage of this approach is however that it introduces an 
ambiguity to the theory, in fact the proportionality coefficient relating mass 
to charge cannot be calculated. Also the very concept of chiral bosonization 
is not very well understood and in passing from fermionic to bosonic for- 
mulation some hand waving arguments are needed (see for example the nice 



papers of Harada |28]). Below, we quantize the chiral Schwinger model in its 
original fermionic formulation and show that the resulting theory is solvable 
and unambiguous. 

The action describing the chiral Schwinger model on the space time of 
topology of X reads 



i_ J d'a:{-^F,,F^'' + + eA,]{l - 7^)^), (6.1) 



where the fiat Minkowski metric rj^u = (— 1,+1) and 7' 







7^=^^ o)'^^~(^ 1)'"'"'^ terms of vector potential components 

Afj^, fi = 0,1 and the chiral component of spinor field x the above action can 
be rewritten as follows 

/ = ^ld'xi^idoA,-d,Aoy-txHdoX + dix) 
- ex^x{A + Ai)) 

= Lem + Lp + Lint- (6.2) 
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The form of Lagrangian suggests that the theory will be greatly simplified if 
we use the light cone coordinates x"^ = ± x^. Indeed in this coordinates 
the Lagrangian reads simply 

Hx^) = ^ / dx~i^id^A+ - d+A^y - ix^d+x - ex\A+) (6.3) 
We immediately see that there are constraints in the problem, to wit 

=0. (6.4) 



d{d+A+) 

Using the expression for momentum of A. 



BT 1 

we can easily calculate the canonical Hamiltonian 

Hc=^j dx-i}-P-p- - A+{d.p- - ex^x)), (6-6) 

which is indeed very simple. After quantizing in coordinate representation, 
P = —i-^ the expression ( |6.6| ) becomes a quantum Hamiltonian operator. 
In the formula above we solved already the fermionic second-class constraints 
i.e., the fundamental fermionic bracket is 

[xHx')^x{y')]+ = -y'), (6.7) 

as usual. From the consistency condition 

[P+,H]=0, (6.8) 

we obtain the secondary constraints 

g{x-) = d^P~ -ex^X = 0, (6.9) 

which is nothing but the Gauss law. 

To proceed let us expand all fields present in the problem into Fourier 
modes: 

A_ = ^One-*""^", 4 = a_„,n>0, 

n 

P- = E^ne-""^ 4 = 7r_„,n>0, 

n 

A+ = ^a^e™^", ai = a.r.n>0. (6.10) 
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We will expand the fermionic field in an unusual way, which is however 
essential for later purposes 



X(x-) = ^ w„e^-exp(-e^ -(a^e-^"^^- - aLe^'"^")), (6.11) 

— Tfl 

n m 

and 

X\x-) = 5:4e-^"^exp(eE -(a„e-"^^" - a^e^™--)), (6.12) 

— — Tfl 

n m 

The use of such expansion can be explained as follows. From purely prag- 
matic point of view, this is the expansion which directly leads to the de- 
sired final identification of the spectrum of the chiral Schwinger model. On 
the other hand, this kind of expansion is to be expected from the following 
physical reasoning. It is well-known that because there are no propagating 
photons in two dimensions, the states of the model should consist of two 
fermions linked by the string (in the sense of the old string model of QCD). 
Then the quanta u and could be interpreted as a fermion with 'half of 
the string'. Finally, it should be mention that the expansion above leads 



to the so-called 'kinematic normal ordering' advocated in p9[ and success 



fully employed to compute anomalous commutators in range of two- and 



four- dimensional theories (see [0 and references therein). 

From these expansions we see that the canonical anticommutational re- 
lations for the fermionic modes holds: 

[Un,ul^]+ = 6m,n- (6.13) 

Also we take 

[nn,al^] = [Kl,am] = -iSn,m- (6.14) 
From the fact that [P_ , x] = [P- > X^] = we find that 

[T^ln^Un] = -e—Un+m, (6.15) 

m 

[vTmjMn] = e—Un-m, (6.16) 

m 

[t^LuI] = e—ui_m, (6.17) 
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[TTm,ui] = -e—ui+rn- (6.18) 

TTli 



In terms of Fourier modes the Hamiltonian (|6.6|) reads 



H = E(^n^l + ^nGn + OWn) + ^4- (6-19) 
n>0 ^ 

Here Qn and Qj^ are Fourier components of the Gauss law constraint which 
are aheady in the Gupta-Bleuler form 

Gi = -^nTxi + e^u\_^uu > (6.20) 

Qn = iuTTn + e'^uj^^ui, n>0 (6.21) 

Q = 5]uS + go. (6.22) 

In the last equation we introduce the normalization constant go for later 
purposes. Let us observe that the constraints Q, commute with Un and 
u\^, which was a reason for introducing the unusual expansion for x and 
above. 

To proceed, us define the fermionic Fock space as follows. First we split 
the operators ui, u\ with respect to the vacuum state |0 >f, 

= rft n < (^-23) 



and 

6„|0 >F= d„|0 >F= 0. (6.25) 

It is easy to see that the charge constraint Q is first-class even on quantum 
level as it commutes with both Qn and Q}^. The algebra of other constraints 
is of the Gupta-Bleuler type, 

[Qn.Qm] = {QlQl]=Q, (6.26) 

[Qn, Ql] = e^n5m,n. (6.27) 
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These formulae are derived in appendix. In the derivation we use operators 
4>n = J2i 'A_nUi and 0jj = i^l+n'^i which are nothing but (anomalous) 
currents and which together with tt operators satisfy the algebra 



(1)1] = n6rr,m, [0n, 0m] = [^l ^L] = 0, (6.28) 

kn,0|„] = [4,0™]=O, (6.29) 

[7^n,4>m] = [T^iAlnl = -'^<^^n,m- (6.30) 

We can easily find a representation of these commutational relations and 



d 

d(f)ii 

, d ie ,1 ie^ .„ „„, 

TTn = -^—J + —4>l - IT^n 6.32 
dak n 2n 

We can also write down the representation for Q and 

Gn = n—j- + —an (6.34) 
oah ^ 

Due to the anomaly the Gauss-law constraints are second-class. Taking Qn to 
be a holomorphic one, which together with Q describes the physical states, 
and proceeding as in Section 4, we find easily the Gupta-Bleuler Hamiltonian 

— 4>n'^n + -ttq + constraints, (6.36) 
n>o 2 

which commutes with Qn by construction. Since ^„ = vr^ + — 0jj the form of 
the Hamiltonian can be further simplified to give 

HgB = -Y1 -^nfl + -VTo' + ^r^Gn + loQ- (6.37) 
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Let us observe that the Hamiltonian (^) and the holomorphic constraints 
Qn are almost identical to the corresponding ones for the scalar field dis- 
cussed before. The only difference is the ItTq term in the Schwinger model 
Hamiltonian and the presence of Q constraint. 

Having constructed all the relevant operators let us turn to revealing the 
physical Hilbert space and the dynamics of the theory. The Q constraint can 
be solved easily by observing that it commutes with and 0^ and annihilates 
the fermionic vacuum state once the constant go is properly adjusted. Then 
the physical states annihilated by Q are of the form^ 

^ = EV^n..m■.M(«,«^«o)(0^)^ ■ (0^)Z:/|O >F, (6.38) 

where the range of summation is finite. Let us now act with the holomorphic 
constraints on the physical states. Since Qn commute with fermionic opera- 
tors and annihilates the vacuum by construction, the condition Qn\phys >= 
reduces to 

d ^ 

(ri— p + jan)iJ = 0, (6.39) 

from which it follows that 

^ = e"^S">o n""""?/;o(a, ao). (6.40) 
Acting on physical states, the Hamiltonian (^ takes a form 

H = e'T. + U. (6.41) 

„>o n d(f)k 2 

The first term in this operator is nothing but the Hamiltonian operator for 
a massive scalar field with mass = i"- To discuss the second term let 
us observe first that the configuration space of is [0, 1] with endpoints 
identified. Indeed there is a remnant of the gauge symmetries of the initial 
theory still present, given by Sao = g~^dg, g = e**^^, G Z, by virtue of which 
one can transform any into above mentioned interval. Therefore the Oq 



^^The states presented below are certainly annihilated by Q. To show that the opposite 
is also true, i.e., that all the states with this property are of the form of (6.41) is more 
difficult. This type of questions has been extensively analysed in the context of string 
theory. A reader may consult e.g. BQ] for details. 
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part of the wave function is 6"^™° and the action of ttq merely shifts all 
the energy levels by the same amount and therefore is physically irrelevant. 
Therefore we conclude that the quantum theory of Chiral Schwinger Model 
with a charge e is equivalent to the family of the theories of massive scalar 
field, with mass = y, parametrized by the eigenvalues of the ttq operator. 

7 The Antianomaly 

The anomaly, in the language of canonical quantization is the appearance of 
an additional central term in the algebra of classically first-class constraints. 
It may happen therefore that for the classical Poisson bracket algebra of 
symmetry generators 

{gi,gj} = -fj/gK, (7.1) 

the corresponding commutator algebra reads 

[GI,GJ] = ^fI/ Gk + uju, (7.2) 

where u is the anomaly. This means that the system which possessed a 
set of local symmetries on classical level, is not gauge invariant quantum 
mechanically. 

The lesson we learned from the sects. 3 and 4 is that using the Gupta- 
Bleuler procedure any theory with a general constraints commutator algebra, 
(|7.2|) can be consistently quantized. Even more, since the procedure described 
above clearly gives the standard answer for systems which can be quantized 
according to Dirac and can be applied also in more general situations, we 
will consider it as our quantization tool. 

Let us now return to the systems with anomalies. We showed above that 
due to quantum corrections some of the classical gauge symmetries may be 
not present in quantum theory. However the opposite situation may occur 
as well: given the classical algebra of the form 

{gi,gj} = -fi/gK + u;u, (r.3) 

the quantum counterpart may be closed, to wit 

[GI,QJ]=^fI/GK. (7.4) 
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I will call this effect antianomaly. From the physical point of view an- 
tianomaly means that the quantum system has larger set of symmetries then 
its classical counterpart. It is even possible that we have a classical theory 
with no symmetries at all, which is becoming a gauge theory after quantiza- 
tion! The example of such theory is presented below. 

Let us remark at this point that no standard technique of quantization 
could be applied in this all of them require the recognition of the 

constraints algebra of the theory already at the classical level i.e., prior to 
quantization. Indeed, when one attempts to quantize a theory described by 
the set of constraints, which form the second class algebra according to the 
Dirac prescription, one uses the Dirac bracket and removes the constraints 
Qi from the theory. There is therefore no way of seeing the effective gauge 
symmetry of quantum theory, described by the algebra (|7.4]) . One may also 
observe that if the number of constraints Qi is N, and the dimension of the 
phase space of the problem is 2D (in the infinite-dimensional case one should 
count per space-point), the theory resulting from the Dirac quantization of 
( |7.3| ) will describe a system with D — |iV physical degrees of freedom, as 
compared with D — N degrees of freedom described by (|7.4| ). These two 
theories are therefore inequivalent. At this point one may argue that we are 
free to choose any prescription, however the explicit model discussed below 
shows that only the second approach is correct. Let us stress two points. 
Firstly, the very discussion of the antianomaly problem is only possible be- 
cause, having in hand the Gupta-Bleuler procedure described above, we are 
no longer afraid to quantize a theory with second-class constraints unsolved. 
Secondly, the problems of the Dirac quantization propagates to the path in- 
tegral formulation of the theory with second-class constraints, which is clear 
from the derivation of such path integral 0]. 

The most important message from the derivation above is that in general 
there is no direct relation between the symmetries of quantum theory and its 
classical counterpart. Therefore one should not use any information concern- 
ing classical symmetries in an analysis of quantum systems, before proving 
that they are identical to the quantum gauge invariances. 

It is not hard to find a physical system with antianomaly. Take the 
(non - chiral) Schwinger model i.e., the two-dimensional QED coupled to a 
single (for simplicity) Dirac fermion. It is well known that this theory is 
anomaly free. Now let us integrate out (in the path integral formalism) a 
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left handed fermion. The resulting effective action will be non - local, but 
the locality may be reestablished by introducing an additional (chiral) scalar 
field into the theoryQ |2^. The theory formulated in this way is, of course, 
completely equivalent to the initial Schwinger model (up to the ambiguity 
mentioned in the footnote) and therefore must be anomaly free. The virtue 
of introducing the scalar to replace the chiral fermion (which is nothing but 
the bosonization of the theory) is, however, that the quantum mechanical 
anomaly has been translated into a classical second-class constraint and there 
is no gauge symmetry on the classical level. Since, as said before, quantum 
mechanically the model has a U{1) gauge symmetry the antianomaly effect 
must occur. 

Let us present the Lagrangian of the Schwinger model with partially 
bosonized fermion field, in the light cone coordinates. 

L{x+) = J dx-{^id+A_ - d^A+f + + 2eK0FA+4v« 

+ i(9+0)(9_0)-^(9_0)2 + ei(a+-9_)0A_ 
- ^elA'_ + elaA+A^}, (7.5) 

where a is the above mentioned ambiguity parameter and to be as general as 
possible we denoted right (left) charges by cr (cl) respectively. It is easy to 
see that there is no gauge invariance in the theory defined by the Lagrangian 
( |7. 5| ) . We will show however that the corresponding quantum theory has a 
local U{1) symmetry. 

It is straightforward to calculate momenta 11$ = ^^^^^^ and then to con- 
struct the canonical Hamiltonian 

n = Jdx-{^P^-A+d^P_ + ^{d^<pf~2eRy/^A+ijyR 

+ eLd^M- + lelAl-elaA+A.}, (7.6) 

There are several technical problems which must be solved in a course of this pro- 
cedure. Not least it happens that the additional parameter appears, which reflect the 
ambiguity in calculation of the functional determinant. This parameter can be however 
fixed by demanding the gauge invariance of the resulting theory, as will be shown later. 
Also an additional outside input is required if one wants to obtain the chiral scalar instead 
of scalar field, see ||2^ . 
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and the constraints 

TTv^fl + #R = 0, (7.7) 

= 114,- ^d^(P - clA. = (7.8) 

P+ = 0, (7.9) 

g = -9_P_ - elaA^ - 2eB^V^i)\il)R = 0. (7.10) 

In formulae above we denoted momenta associated with the electromagnetic 
field by P^. The first constraint ( |7.7| ) is the standard second-class con- 
straint for fermions, reflecting the fact that the kinetic term for these fields 
is linear in the time derivative. Since it does not change in course of quan- 
tization we simply solve it by using the Dirac bracket. Then we find the 
fundamental anticommutator 

[i;^j,ix),My)]+ = Six-y). (7.11) 

The constraint ( [7.9|) is clearly first-class and refiects the fact that the wave 
function should be A^-independent. The remaining two constraints, if treated 
classically form a second-class system, to wit 

{g{x),g{y)} = -2elaS\x-y) (7.12) 
{C{x),Ciy)} = -5\x-y) (7.13) 
{g{x),C{y)} = -eL6'{x-y) (7.14) 

Before proceeding let us recast the constraints into simpler form. First, by 
making the canonical transformation vr^ — > vr^ — clA^^P^ — P_ — 6^0 we 
turn the constraint (|7.8|) into the standard chiral boson constraint 



= 714- U^<p. (7.15) 

Then we diagonalize the constraints system by taking a linear combination 
of them, 

g ^ g - eiC 

= -9_P_-eL(7r^ + ^9_0)-eiaA_-2eKv^^Jj^i?, (7.16) 
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such that the algebra of new constraints reads 



{<3{x)My)] = el{l-2a)5'{x-y), (7.17) 
{C(x),C(y)} = (7.18) 
{Q{x),C{y)} = 0. (7.19) 

(7.20) 

The dynamics of the theory is governed by the primary Hamiltonian 

Hp = |rfa;-{^(P_+ei0)2 + l(9_</.)2 + ei9_</.A_ + leiA! 

+ A+g + \C}, (7.21) 

where A is an Lagrange multipher. In the evaluation of these formulae we 
took into account the fact that the Poisson bracket of fermionic bilinears 
ipR'4'R is equal to zero. This, however, changes in the quantum theory - 
indeed it is well known that the commutator of fermion bilinears is 

[4V^«(x),4^«(l/)] = - y) (7.22) 

Taking this into account we see that the quantum commutator of the 'Gauss 
law' constraint Q reads 

[g{x), g{y)] = z6'{x - y)[el{l - 2a) + el]. (7.23) 

We can also calculate the commutator of the constraint Q with the primary 
Hamiltonian, 



[g{x),np] = J dy{tei{a - 1){P^ + eL<j)){x)S{x ~ y) 

+ tA^S'ix-y)[elil-2a) + el]}. (7.24) 

We see therefore that the 'Gauss law' constraint can be made first-class and 
therefore a generator of gauge transformation if a = 1 and e| = ej^. For 
these values of parameters the model described by the Lagrangian ( [7.5|) is 
equivalent to the standard Schwinger model and, as we just have shown, 
exhibits the antianomaly phenomenon. 

We discovered above that the parameter a should take the value 1. There- 
fore the theory described by Lagrangian ( |7.5| ) with a = 1 and with no right 



77 



fermions should be equivalent to the fermionic theory discussed in Section 
6. To check this let us look at the equation of motion following from this 
Lagrangian. 

Assuming appropriate behaviour of all fields at spatial infinity, from the 
field equation (we drop the subscript L) 

(9+ -9_)(9_0 + eA_) = 0, (7.25) 

we learn that 

a„0=-eA_. (7.26) 

Substituting this result into the A_ equation of motion we see that the com- 
bination of fields p = ^9+0 + is a massive scalar field: 

d+d^p + eV = 0, (7.27) 

with mass = e^. Taking these two results into account the equation for 
reads 

(9!p + ea_0 = O. (7.28) 

Acting on the right-hand-side of this equation by 9+ and using (|7.27| ) we see 
that 

9_(ep-(9+0) = 9_A+ = 0, (7.29) 

which is the equation defining a chiral scalar field. Therefore the spectrum of 
the left chiral part the model consists of one massive and one massless chiral 
scalar field and agrees with the spectrum of the chiral Schwinger model with 
left handed fermions described above. 

Let me conclude this section with some remarks. The model described 
and solved above is of course artificial and its sole role was to give an example 
of the antianomaly phenomenon. However, as an extra bonus it provided us 
with a way of deriving the actual value of a parameter. It remains still an 
open question whether there exist less trivial (e.g. four dimensional) theories 
with antianomalies. The constraints structure of such theories is of course 
pretty easy to guess: for example one can take ^ in (11) to include a fermionic 
chiral current and then to arrange uu to be —i times the chiral anomaly. 
The hard part however would be to construct an appropriate Hamiltonian, 
leading to unitary, Lorentz invariant and hopefully renormalizable quantum 
theory. Finally, let me stress that antianomalies have nothing to do with 
anomaly cancellation, which is an arranging for such a field content of the 
theory that the classical and quantum symmetry groups are identical. 



78 



8 Conclusions 



Let me briefly summarize the contents of this paper. In the first three sec- 
tions I derived the canonical Gupta - Bleuler approach to quantizing of the 
anomalous theories. I presented a 'cook book' describing how one should 
proceed to quantize such a theory from the scratch. In Section 6 I used this 
procedure to quantize the 2 - dimensional Schwinger model and in Section 7 
(among other things) I showed that the result coincides with the well- known 
standard one obtained in the past. 

In my opinion, however, the most important result of this paper is the 
derivation of additional terms present in the Hamiltonian path integral for- 
mulations of anomalous theories. It may happen that the path integral in 
the form presented in this section can be applied to consistent quantization 
of anomalous theories in the physical dimension d = 4 like chiral QED or 
QCD. It would be interesting to see what the predictions of such theories 
are, for example if the anomalies give rise to mass generation, as they do in 
two dimensions, and if some realistic models describing elementary particle 
physics can be generated in the framework of such theories. Work in this 
direction is in progress. 
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9 Appendix 

In this appendix we collect some results of more technical nature which are 
used in the main text. 
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To start, let us denote for integer n (pn = Hiu\_nUi and = J^i^jj^n'^i. 
Let us now calculate the commutator algebra of these objects: 



k,l 
/ 

The sum in the last line of above expression is not well defined and requires 
regularization. We do that by taking the range of summation in definition 
of (pn and 0jj to be finite and range from —L to L where L is large enough. 
Using this regularization, the last line above becomes 

L 

Y,{u]_^Ui_m - uj^^Ui+n). (A.l) 
-L 

li n ^ m one can shift the summation range in the second term to see that 
all terms can cancel when L goes to infinity and the sum is zero. For n = m 
however we have products of non - anticommuting operators and more care 
is needed. Acting on states with finite number of fermions, such that only 
states with b}^, d\^, L > n,m are present, we rearrange both sums to the 
normal ordered form and see therefore that 

L 

'^{ul_^Ul-ra - ul_^^Ul+n) = n6n,m- (A.2) 
-L 

We stress here that the appearance of the non - trivial right-hand-side in 



(|A.2|) is a purely quantum effect related to the ordering problem, as opposite 
to the classical theory (in particular the Poisson bracket of (j) and (p^ vanishes). 
It should be also stressed that this formula gives unambiguous result only 
on the Hilbert space consisting of finite number of fermionic excitations, for 
infinite number of fermions it again becomes ambiguous. However, we extend 
our result also to the wave functions with infinite number of fermionic states, 
which will be needed in description of physical states. Actually the non - 
vanishing of the RHS is responsible for the anomalous term in the algebra of 
classically first-class constraints Qn and Qj^, as will be seen below. 
In analogous fashion one can prove that 

[0„,0m] = [0t,0]n] = O, (A.3) 
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K, 0m] = [vr„, (pU = 0, 



(A.4) 



kn) 0m] = K, 0m] = "^^n.m- 

Now we are in position to calculate the algebra of constraints: 



(A.6) 



(A.5) 



It is easy to see however that 



[Gn, GU = 



= ne^Sn, m. 



(A.7) 



which is the anomaly of the Gauss law constraint. 
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